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Abstract 

We have investigated the nonlinear amplitude vector equation gov- 
erning the evolution of optical pulses in optical and UV region. We are 
normalizing this equation for the cases of different and equal transverse 
and longitudinal size of optical pulses or modulated optical waves, of 
week and strong dispersion. This gives us the possibility to reduce 
the amplitude equation to different nonlinear evolution equations in 
the partial cases. For some of these nonlinear equations exact vortex 
solutions are found. Conditions for experimental observations of these 
vortices are determined. We have found large spectral zones on the 
'wings' of the electron resonances of metal vapors and also near to 
the plasma frequency where the amplitude equation can be reduced 
to 3D+1 vector nonlinear Schrodinger equation. 

PACS 42.81.Dp;05.45.Yv;42.65.Tg 



1 



1 Introduction 



The investigation of optical pulses and of modulated periodical waves in wave 
guide structures as fibers and planar waveguides is connected with the addi- 
tional requirements to the optical property of the materials. One standard 
way for confinement the optical beams and pulses is by variation of the linear 
or nonlinear refractive index in plane, which is transverse to the propagation 
direction. In an isotropic dispersive media the refractive index depends only 
on the frequency, the group velocity depends on the chromatic dispersion and 
the wave guide modes can not be introduced in linear regime of propagation. 
To provide one correct analysis in this case we need to get in mind also the 
different spatial size of the optical pulses. The longitudinal spatial size is de- 
termined from the time duration by the relations zq = v g to = cto/n g , where 
a typical value of the group index n g in transparency region of pure silica is 
{rig = 1.5). For example, the typical transverse size of laser pulses is evalu- 
ated from rj_ = 1 — 5 mm to r_|_ = 100 fim. While for the longitudinal size of 
pulses or spatial period of modulated waves, there are mainly two possibili- 
ties: a. Pulses from nanoseconds up to 40 — 100 picoseconds. For such pulses 
zq » r_i_ and their shape is more close to optical filaments, b. Pulses from 
few picoseconds to 100 femtoseconds. In this case zo ~ r± and the pulses look 
like as optical bullets. When the transverse and longitudinal part are approx- 
imately equal (case b), the ignoring of the second derivative in z direction in 
amplitude equation is not possible. The propagation of such type of short 
optical pulses in nonlinear media must be characterized, not only with the 
equal dimensions in x, y and z direction, but also with non-stationary optical 
linear and nonlinear response of the media. That is why we can expect, that 
the equations governing the propagation of such type of pulses, like optical 
bullets, are different from the well known paraxial approximated equations, 
which describe evolution of nanosecond and picosecond pulses [TJ |21 El • An- 
other basic characteristic of the pulses is their polarization. As it is well 
known, in the case of plane monochromatic wave, the electromagnetic field is 
polarized in the plane perpendicular to the direction of propagation. In this 
case four Stokes parameters, two amplitudes and two phases, connected with 
the transverse components of the electromagnetic field can be introduced 0] . 
This corresponds to SU(2) symmetry of propagation of the monochromatic 
field. The nanosecond and hundred picosecond optical pulses admit very 
little spectral bandwidth, and the polarization component in the propaga- 
tion direction is small in respect to the transverse components. With good 
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approximation these pulses or modulated waves also can be investigated as 
polarized in x, y plane. Their spatial form also is more close to the cw wave, 
as the longitudinal size is hundred times greater, than the transverse size. 
In contrast to the case of long pulses, the few ps and femtosecond optical 
waves are with relatively equal size in x, y, z direction and admit large spec- 
tral width. For them such coherent SU(2) tensor can not be introduced. The 
experimental results show, that the electric field of sub-picosecond and fem- 
tosecond pulses, always contains also third longitudinal component. With 
the enlargement of the spectral width of the vector field also the longitudinal 
component grows (the component normal to the standard Stokes coherent 
polarization plane). Recently Carozzi et al. have shown that for the wave 
packets higher class of symmetry - SU(3) exists [5j. Using the presentation 
of Gell-Mann of SU(3), they prove that five independent parameters, three 
amplitudes and two phases, define the dynamics of propagation of the vec- 
tor field. This corresponds to a three-component vector field. The above 
made analysis and the experimental results show that in the investigation 
of sub-picosecond and femtosecond pulses we must always have in mind the 
three-component vector character of the electromagnetic field. An initial in- 
vestigation of the propagation of the optical pulses was provided in the frame 
of spatiotemporal scalar evolution equation [T3 [THI - We found that vector 
generalization of such spatiotemporal model can be applied to optical pulses 
or to modulated periodical waves in UV region for dielectrics with strong 
dispersion. 



2 Nonlinear amplitude vector equation 

Starting from the Maxwell equation for an isotropic, dispersive, nonlinear 
Kerr-type media with no stationary response, we derive the next slowly vary- 
ing vector amplitude equation [TE] : 
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where k , v, k' Q and n 2 are the carrying wave number, group velocity, dis- 
persion of the group velocity, and the nonlinear refractive index respectively. 

— # 

With A we denote the slowly varying vector amplitude of the electrical field. 
The equation ((TJ is written in second approximation to the linear dispersion 
and in first approximation to the nonlinear dispersion. The kind of nonlinear- 
ity is connected with the fact, that in this paper, as in the previous one [To] , 
we investigate only linearly polarized electrical field. The dependance of the 
nonlinear polarization from different kinds of polarizations of the electrical 
field is discussed in the Appendix 1. We will apply to the basic amplitude 
equation (JTJ) the method of different and equal transverse and longitudinal 
spatial scales, as it is connected with the natural initial shape of the optical 
pulses. This application is one simplification of the multi-scale method, in- 
troduced in the nonlinear optics in [SJ Q |S] and also in the hydrodynamics 

m 

3 Case of different transverse and longitudi- 
nal size 

Applying a "moving in time" (t' = t — z/v;z' = z) transformation to the 
vector amplitude equation (JTJ), we obtain: 
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where Aj_ = ^ + -§^p- The estimation of the influence of different terms 
in equation (0) can be reached by writing it in dimensionless variables. To 
make one assessment of the difference between transverse and longitudinal 
dimension of the pulses we are introducing separately transverse r± and lon- 
gitudinal Zq constants. Defining the rescaled variables: 



A = A A'; x = r±x'; y = r±y'; z' = z Q z" ; t' = t t" . 



(3) 
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and constants: 



W*o = v; a = k z ; 5 = r ± /z ; z dis = t /k"; 



Pi = k r ± / Zdi S ; 7 = k r ± n 2 \Ao\ /2; 71 = n 2 |A | /2, 



(4) 



equation (J2J) can be transformed in the following (the primes and the seconds 
are not written for clarity): 
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where = t^/k" and in the transparency region we use the approximation 
dn 2 /du ~ 0. 



3.1 Pulse propagation in optical region: week disper- 
sion 

In the optical region the wave vector is valued from ko ~ 10 4 to k^ ~ 10 5 
cm -1 . The typical value of the dispersion parameter in the transparency 
optical region of dielectrics is from fc" ~ 10~ 28 s 2 /cm up to k" ~ 10~ 24 s 2 /cm 
in UV region |17j . The transfer sizes r_i_ of laser generated optical pulses or 
modulated periodical in time waves are valued from r± < 1 to r± ~ 10 -2 
cm. We will investigate laser pulses or modulated periodical waves with 
time duration or modulation period from few nanoseconds (to ~ 10~ 9 s) up 
to 40 — 100 picoseconds (t ~ 4 — 10 x 10~ n s) in the transparency region of 
a nonlinear Kerr type media. The group velocity index in this case is about 
n g ~ 1.5 and using the relations z = v g t = ct /n g ~ 10 1 — 10 2 cm, we find 
that the dimensionless parameter in front of second derivative in z direction 
and cross-term is very small 5 2 = r\jz\ ~ 10~ 2 — 10~ 4 . Another important 
relation is that if slowly-varying amplitudes are used a = k zo ~ 10 3 — 10 4 . 
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For such typical values of the dispersion k n , wave vector fc , time duration i 
and transverse dimension r_|_ of the optical wave, the dimensionless parameter 
Pi in front of the dispersion term in equation (JSJ is very small P = for\j Zdi S ~ 
10 -3 — 10~ 4 . Here we study the cases with power near the critical for self- 
focusing 7 = 1. These valuations give us possibility to estimate dimensionless 
constant in front of the different terms of normalized amplitude equation (0) 
for the transparency optical region of a dispersive nonlinear Kerr type media: 

a = k z ~ 10 3 - 10 4 ; S 2 = r\/ z 2 ~ 1(T 2 - KT 4 ; 

Pi = fori/ Zdis ^ 10~ 3 - 10" 4 ; (6) 
7 = k 2 r 2 ± n 2 \A \ 2 ~ 1; 7l = n 2 \A \ 2 < 1. 

Neglecting the small terms in (J3J) we obtain the next amplitude equation: 



f)A i _ -> 
-i— + -A ± A + A A = 0. (7) 

oz 2 

As it can be expected, the dynamics of such kind of long optical pulses 
is governed by the well known paraxial approximation. This equation was 
introduced initially for optical beam PHI21 El- The fact, that with such type 
of nonlinear paraxial equation it is possible to investigate non-stationary 
processes is well known also fUl El Ej ■ It seems natural the dynamics of 
propagation of long pulses to be governed by the equation which describes 
dynamics of optical beam. The paraxial approximation is applicable even 
when the relation between the longitudinal and the transverse part is of the 
order of ten (5 = r±/z ~ 1/10). The dimensionless coefficient in front of 
the term with second derivative in z direction and cross-term is of the order 
of square of 5 ( 5 2 ~ 1/100 << 1), the term can be neglected and this 
leads again to the paraxial equation (JJJ). Effects of self-steepening due to 
the first order of the nonlinear dispersion (second term in equation (j3J)) is 
not possible as always 71 ~ n 2 \ A | 2 << 1. Different kind of generalization of 
equation ((7j) in respect to nonlinearities [121 HH1 UH 120] and multi-component 
vector fields |2*T1 122] was performed. It is important to point here that the 
second derivative in z direction and the crossed term are neglected from the 
difference between transverse and longitudinal size (5 2 = r\j 'z\ ~ 10~ 2 — 
10" 4 << 1), and not because the slowly varying amplitude is used. When 
we investigate the optical pulses with equal transverse and longitudinal size 
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5 2 = t\/zq ~ 1, the second derivative in z direction and cross term of the 
amplitude function are of the same order that the transverse Laplassian, and 
these terms influence considerably on the dynamics of propagation. This 
important case we will investigate in the next section. When the power of 
the pulse is less than the critical for self-focusing 7 << 1 the diffraction term 
dominates and the propagation of the optical pulses is governed by the next 
linear paraxial equation: 

dA 1 A T 

3.2 Optical pulses in UV region: strong dispersion. 
Vortex solutions 

In UV region of media with high density as liquids and dielectrics the disper- 
sion k" increases considerably and can reach the values from fc" ~ 10~ 24 — 
10 -25 s 2 1 cm. The carrying wavenumber is valued from k ~ 10 5 cm" 1 to 
k Q ~ 10 6 cm" 1 and in this spectral region we can reach value with equal 
diffraction and dispersion length, which corresponds to (3\ = k^r\j ' z^ s ~ 1. 
We use for example bulk fused silica and calculate the linear refractive index 
from the Sellmeier relations [2*3] : 

n 2 = 1 + y2 J¥^ ; (9) 

with coefficients B x = 0.6961663, B 2 = 0.4079426, B 3 = 0.897479, wave- 
lengths Ai = 0.0684043 fim \ 2 = 0.1162414 fim A 3 = 9.896161 /im where 
Xj = 2irc/ujj. The choice of carrying wavelength for the optical packet at 
A = 0.264 /im, (w = 0.7131606.10 16 Hz) is connected with the facts, that 
fused silica is still transparent on this wavelength and it can be reached easy 
way by fourth harmonics from Nd YAG laser on A = 1.064 \xm. Considering 
the case of pulses or modulated waves with intensity near the critical for 
self- focusing 7 ~ 1, time duration t — 30 — 40 ps and using the relations 
k{u) = y t uj 2 n 2 '(td) I 'c 2 , fc" = jj^s, the next values for laser characteristics and 
material constants are obtained: 

k = 2.378 x 10 5 cm- 1 ; k" = 1.99 x 10" 25 s 2 /cm; n g ~ 1.65; 
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r ± ~ (1.4 - 1.6) x 10" 1 cm; z dis ~ (0.45 - 0.8) x 10 4 cm; (10) 
z ~ (6 - 8) x 10" 1 cm; fcgrl ~ (1.1 - 1.5) x 10 9 ; 

n 2 |A | 2 ~ (1.1 - 1.5) x 10~ 9 . 

These parameters correspond to the next dimensionless constant: 



A = Vl/^s 10 4 ; 
<5 2 = r^ 2 ~ 10- 2 ; 7 = fegr^na \A \ 2 ~ 1; (11) 
a «5 2 7l = a 5 2 n 2 |A | 2 ~ 10~ 6 < 1. 

Neglecting the small terms from the amplitude equation (JSJ) we obtain the 
next equation governing the propagation of optical pulses with different trans- 
fer and longitudinal size in UV region of dielectrics: 



. r2 dA 1 -> 1 d^A t 2 -. n 

87 + 2 M "2% +1 H (12) 

The scalar variant of this spatiotemporal equation ()12j) was introduced for 
first time by Silberberg ^2] and attracts a growing interest. In the beginning 
this equation was suggested to discover the dynamics of so-called "light bul- 
lets" (LB), i.e. optical pulses with relatively equal transfer and longitudinal 
size. Our investigation corrects this idea and finds that the equation (|T2*|) 
governed the propagation of optical pulses with longitudinal size hundred 
times greater than the transfer one in UV region of dielectrics and liquids. 
In the case of optical waves with very small intensity the nonlinear term can 
be neglected 7 << 1, and we obtain the linear version of the equation (|12j) : 



. r2 dA 1 -j> 1 d 2 A 

Let investigate the two component vector field polarizable in transverse x, y 
plane. For long pulses, as they are more close to the cw regime, the longi- 
tudinal vector component is small and such approximation is possible. We 
represent the optical vector field by two orthogonal components A\x + A 2 y 
and the equation ()12|) is transformed to the next nonlinear system of equa- 
tions: 
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(14) 

-rf^ 2 - + U ±A2 _ 1 A ^ + 7 (L4 2 | 2 + |A| 2 ) A 2 = 0. 



In the linear case the (fT3j) i s reduced to two equal scalar equations for A\ 
and A 2 : 



■ ?2 dA i 1 . . 1 a <9 2 A „ 

where z = 1,2. We rewrite the linear scalar equation (JTKJ) in cylindrical 
coordinates: 



- ia<y 97 + 2 (r ^7 + ~w + ) ~ 2^~dF = °' (16) 

where r = y/x 2 + y 2 and = arctan(x/y). The linear equation (fT6|) admits 
infinite number of exact analytical solutions of kind of: 



At = r n [a„ cos(n6>) + a n sra(n0)] exp (z(/3iz + ^cnR)) , (17) 

= r~ n [a„ cos(n6>) + a n sin(n0)] exp (z'(/3iz + V2a5 2 t)) , (18) 

where n = 0, l,2...oo is one infinite number. The solutions (f 17|) . (fT%|) admit 
modulation frequency f2 = \j2ab 2 and modulation wave vector = /3i in z 
direction. The existence of such kind solutions for monochromatic light was 
observed for first time by Nye and Berry . The solutions of kind (JT7jl , (fT5j) 
admit amplitude and phase singularities. One elegant way to remove the 
amplitude singularities by using the combination of Gaussian envelope and 
the solution (|18j) with n — 1 was suggested for first time in for paraxial 
beam. This method is applicable also to the linear equation (|16|) and optical 
vortices in combined pulses with finite amplitudes can be found. Usually the 
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localized solutions of the linear equations are not stable. Now we will turn 
to most interesting case; exact solutions of the nonlinear system of equation 
([14)1 . Again we rewrite the system (|T4")) in cylindrical variables: 

_ ia p dAi + 1 ( 19Ai + + 1 92 A A 

dz 2 \r dr dr 2 r 2 d9 2 J 

-\^ + i{\M 2 + \M?)A 1 = ^ 

(19) 

ia5 2dA2 1 / 1 dA 2 (?A2 1 d 2 A 2 \ 
dz 2 \r dr dr 2 r 2 d9 2 J 

The nonlinear system of equations (|19|) admits another infinite number of 
exact vortex solutions of kind of: 

Ay = ^((n + l) 2 - l)/7(r _1 ) cos((n + 1)6) exp (i(f3yZ + V2^6 2 t)) , (20) 
A 2 = ^((71 + 1)2-1)/^-!) S in((n + 1)9) exp (i{(3yz + v 7 ^^) , (21) 

where n = 1,2. ..00 is one infinite number. One important consequence of 
these solutions is, that they are obtained as a balance not only between 
the diffraction and the nonlinearity, but mostly by the nonlinearity and the 
angular distribution. That is why, we can expect one stable propagation of 
such vortices in z direction. Right away we can point here, that conditions 
for finiteness of the energy of the above solutions have been not found up to 
now. 

4 Pulses with equal transverse and longitu- 
dinal size r±_ ~ zq 

The investigation of the propagation of optical pulses or modulated waves 
with equal transverse and longitudinal size r± ~ zq is more convenient to 
provide in coordinate system moving with group velocity. That is why we 
apply a Galilean transformation to the basic vector amplitude equation (£' = 
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t, z' = z — vt), as difference from the case of long pulses, where the "moving 
in time" transformation is more natural. The equation (Q) in this case is 
transformed to: 
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where Aj_ = J^j + We use again the approximation 8712/ du ~ 0. To 



(22) 



estimate the influence of the different terms on the propagation dynamics 
we rewrite the equation ()22j) in dimensionless form. Defining the rescaled 
variables: 



A = A A'; x = r x'; y = r y'; z" = r z'; t" = t t', (23) 
and constants: 



a = 2k r /t v = 2k r ; (3 = v fe" k ; 7 = k r n 2 \A \ ; 



ro/to = v; 71 = 2k r n 2 \A Q \ 2 , 



(24) 



equation (|2*2*|) can be transformed in the following (the primes and the seconds 
are not written for clarity again): 
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We point out, that the normalizing is with equal longitudinal and transverse 
dimensions r . As it was pointed above, it is fulfilled usually for few picosec- 
ond and femtosecond pulses of a pulse with transverse size of 1 — 5 millimeters 
to hundred micrometers. 

4.1 Equation in optical region of media with week dis- 
persion (3 « 1 

Dimensionless dispersion parameter (3 in transparency region of dielectrics 
and gases usually is of the order of 1CT 2 — 1CT 3 . For example we calculate j3 
in the transparency region of pure silica on A = 1, 55/xm using the standard 
Sellmeier expression for dependence of the dielectric constant from frequency, 
and the value really is of this order (3 = — 0, 0069. The constant a has a value 
of a ~ 10 2 (a ~ 2r k ) if the slowly varying approximation is used. In this 
case, when we investigate the optical pulses with power near to the critical 
for self-focusing 7 « 1, 7! < 1, our vector amplitude equation (J25J) can be 
transformed to the next evolution equation: 

The equation ()26|) is one generalization of the well known paraxial equation 
(J2J) in optical region for pulses with equal transverse and longitudinal size. 

4.2 Equation in UV transparency region. Strong dis- 
persion (f3 « 1) 

The critical dimensionless parameter /3, which determines the relation be- 
tween dispersion and diffraction in the normalized amplitude equation ([25)1 . 
as it was established in the previous paragraph, is very little in the optical 
region. In spite of this, the temporal effects are presented in the correspond- 
ing equation and for pulses with equal transverse and longitudinal size 
can not be neglected. The situation in UV region is quite different. In this 
region we can reach (3 ~ 1. When we investigate optical pulses with power, 
near to the critical for self-focusing, the corresponding nonlinear dimension- 
less parameters admit values 7 ~ 1, 71 << 1. Than the equation (J25J) is 
transformed to: 
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dA 7 n d 2 A ln Jd 2 A d 2 A\ T 2 t 

_ !Q _ + A ^- /3 — + -2— j +7 H A = Q. (27) 

If we compare the equations in optical (|2l)Jl and UV region (|27|). it is seen that 
in UV region we can take in mind in addition a negative second derivative 
in z direction of amplitude function and we multiply the temporal effects of 
second order by factor [3 + 1. We use for example again bulk fused silica and 
calculate the linear refractive index n 2 (uj) from the Sellmeier relations. The 
investigation is provided for the spectral region from wavelength A = 0.264 
fim to wavelength A = 0.164 fxm. The group velocity v g , dispersion k" (u) and 
the dimensionless parameter /3(u) are calculated. The typical characteristic 
of the laser pulse are t = 5 x 10 -12 s, group velocity v g ~ 1.8 x 10 10 
cm/s and r± = zq ~ 1 x 10 _1 cm. On Figure ^ Figure El are presented 
correspondingly the graphics of dispersion of the group velocity k" (A) and 
dimensionless parameter /3(A). As it is seen from Fig. 2, the dimensionless 
parameter (3, really can reach values equal to one in UV region for some 
dielectrics. 



4.3 Equation in media with strong negative dispersion 

Now we come to the most interesting case in this investigation: the possibility 
that the temporal effects of second order can be neglected. As it can be traced 
out from basic amplitude equation (0), and to see in the normalized equation 
(l2*5Jl . when k v 2 k" = (3 = —1, the dispersion effects of second order do not 
present in the corresponding amplitude equations and do not influence on the 
propagation dynamics. Recently in ^B] it was found that near the electronic 
resonances of gases and dielectrics and also near the Langmuir frequency in 
cold plasma the dispersion is negative, the dispersion parameter arises rapidly 
and is of the order of (3 m —1. Again we investigate the propagation of optical 
pulses with power near the critical for self- focusing (7 « 1, / y 1 << 1). The 
amplitude equation (J25|) in this case can be written in the next simple form 
of 3D+1 Vector Nonlinear Schrodinger equation (VNSE): 



dA . r d 2 A 



A 



A = 0. (28) 
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Figure 1: Plot of the dispersion of the group velocity k (u) in the UV region 
from A = 264 nm to A = 164 nm for bulk fused silica. 



Here we will study the possibility to obtain (3 = — 1 in the next three cases: 
1. Near the isolated resonance of mercury vapor on A = 0.3653/im. 2. Near 
the Langmuir frequency in cold electron plasma. 3. In deep UV and Ro 
region for dielectrics, gazes and metals. 

4.3.1 Dispersion parameter (3 near the isolated electron resonance 

When e(uj) ~ 1/u the wave vector can be written as: 



k(u) = (29) 



Calculating (3 = v 2 k"k for dispersion relation of kind of ([29)1 . right away 
we obtain (3 = — 1. This gives us a confidence to look for such zones near 
the electronic resonances where e(uj) ~ l/u. Near to some of the isolated 
electron resonances of metal vapors the dielectric constant can be expressed 
by |23| : 
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Figure 2: Plot of the dimensionless parameter (3 = k v 2 k" within wavelengths 
from A = 264 nm to A = 164 nm for bulk fused silica. The dispersion 
parameter reach a value j3 ~ 1 in UV region. 



/47riVe 2 - g 
e(u;) = H ^ . 30 

The investigation are provided for resonance frequency of mercury vapor on 
c^o = 2nc/\ = 0.515568. 10 16 Hz, oscillator strength / = 0.3, and natural 
(i.e. radiative) linewidth 7 = 1/4,53Aq = 1,654.10 s Hz. We use Lorenz 
shape of the line, as we investigate the rarefied vapor and looking in zones 
outlying from resonance on distance more than 47, where the natural shape 
of the line dominated. In Figure 3 it is shown the form of the dielectric 
constant for number density iV = 1, 2.10 15 . In the Figure 4 it is presented 
the dimensionless dispersion parameter (3. As it is well seen on the wing on 
the resonance on distance about 47 we have spectral zone where j3 ~ — 1. 
Right away it must be pointed that (3 depends strongly on the number density 
N. The number density of the order of N ~ 1.10 13 — 1.10 15 appears to be 
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Figure 3: Frequency dependence of the dielectric constant for resonance 
frequency of mercury vapor of u Q = 2nc/\ = 0. 515568. 10 16 Hz (A = 365 
nm), and natural linewidth 7 = 1/4, 53Aq = 1, 654. 10 8 Hz. 

most suitable for large zones with f3 ~ — 1. Seems, for such number density 
the wing of the resonance in best way can be approximated with e(u) ~ 1/u. 

4.3.2 Dispersion parameter (3 near to Langmuir frequency in cold 
plasma 

Dielectric constant in the case of cold electron plasma is given by the expres- 
sion |2"5] : 

e(«) = l-4, (31) 

where uo p = \J AnNie 2 /m is the plasma frequency and N\ is the number of free 
electron particles per cm -3 . Here we investigate a typical number which char- 
acterized the plasma density in generators used for nuclear fusion. In Figure 
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Figure 4: Plot of the dimensionless parameter (3 = k v 2 k" near to resonance 
frequency of mercury vapor on u = 2nc/\ = 0. 515568. 10 16 Hz. On distance 
47 we find the spectral zone where j3 ~ — 1. 



5 is presented s{uj) for frequencies higher and close to plasma frequency and 
number density N\ = 1.10 17 cm -3 . The frequency is normalizing on spectral 
width of optical pulse 5oo = 1.10 10 Hz, which corresponds to pulse with time 
duration At ~ 2.10 -10 seconds. The corresponding to this dielectric constant 
dispersion parameter /3(u) is presented in Figure 6. As it is seen from the 
graphics, the large spectral zone exists, more than twenty spectral widths 
from plasma frequency, where (5 = — 1. The transverse and longitudinal size 
such optical pulse valued from 0.5 — 1 cm. An important relation is pointed 
out with investigation of dependence of /3 from time duration of the pulses 
and number density. With decreasing the time and spatial duration of the 
pulses (fs pulses) and increasing the number density, this large spectral zone 
with (3 = —I keeps on. 
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Figure 5: Frequency dependence of the dielectric constant near to Langmuir 
frequency in cold plasma for number density Ni ~ 1.10 17 . The frequency is 
normalizing on spectral width of the optical pulse 5uj = 1.10 10 Hz. 



4.3.3 Dielectric constant e(u) and dispersion parameter (3 from 
deep UV to R6 region 

As it was established in j2E|, it is possible to find equal dispersion relation 
of the dielectric constant e(u) for different kind of materials, as dielectrics, 
gases, semiconductors, metals, when the carrying frequency of optical wave 
is much higher than the resonances of this media. In these cases the localized 
energy of optical wave is much higher than the critical for ionization of the 
materials and the linear and nonlinear response is the same as in cold plasma. 
Practically the electrons interact as free electrons on these frequencies. It is 
not hard to find that in this region again we have the dispersion relation of 
dielectric constant of the kind of: 



e(u>) 



1 

to 2 ' 



(32) 
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Figure 6: The corresponding dimensionless dispersion parameter (3 near the 
plasma frequency. Large spectral sone with (3 ~ — 1 is established. 



The critical parameter which characterized the plasma frequency in these 
materials is the number density N\ of the free electrons. The typical number 
for semiconductors is N\ = 1.10 18 and for metals is N± = 1.10 22 . That is why, 
the application of this formula started from deep UV for light elements as H, 
Li, e.c. and from R6 frequencies in heavy elements. For optical waves when 
uj > uj p , there are no difference between metals, semi-conductors, dielectrics 
and metals. All of them are transparent for optical waves with carrying 
frequency u > u> p . Again we can obtain, that near the corresponding plasma 
frequency of these materials, large frequency zone with negative dispersion 
parameter of order of (3 ~ —1 exists, if the picosecond and femtosecond 
pulses are used. For all of these three important cases the propagation of 
localized optical pulses is governed by the vector 3D+1 Nonlinear Schrodinger 
Equation ([28)1 . The numerical experiment provided in JH] show a stable 
propagation of the vortex solutions. The intensity picture of these multi- 
component solutions is truly symmetric and look like as stable Fraunhofer 
distribution. The internal structure of these vortices is presented by the 
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spherical harmonics with integer number. 



5 Exact vortex solutions on VNSE 

The solutions of the 3D+1 Vector Nonlinear Schrodinger amplitude equation 
(|2HJ) in a fixed basis are [16] : 



A, 



^exp \ i\J a£lr 
V2 * 



sin 9 cos ip exp (iVtt) 



(33) 



exp (iy/aTlrj 
A y = v 2 sin 9 sin ip exp (ifit) , 



(34) 



A 7 



^exp \i\J aVLi 
V2 * 



cos 9 exp (iQt) 



(35) 



where r = \/x 2 + y 2 + z 2 , 9 = arccos f, and p = arctan^ are the moving 
spherical variables of the independent variables x, y, z. The corresponding 
real solutions ^t(A), are: 



1 ^ e 
U(A X ) = - (A x - A* x ) = V2- 



sin ( VaOr + Qt 



sin 9 cos <p, 



(36) 



HA y ) = j i {A y - a; 



^sin ( y/aQr + Fit 
V2 i ■ 



sin # sin y?, (37) 



1 ^sin ( y/aQr + Clt) 

&(A Z ) = -(A z - A* z ) = V2 ^ L cos 9. 

2t r 



(38) 



It is important to point, that in \%jg^ conditions for finiteness of the energy 
of the solutions (|3*B |l -(|3*5 jl was found. As it can be established, the averaged 
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Figure 7: Fraunhofer distribution for the averaged with one time period 
intensity of the solutions (|39|) of 3D+1 Vector Nonlinear Schrodinger equation 

(123). 



with one time period intensity field of these solutions admits total radial 
symmetry: 



(A x ) 2 + (A y ) 2 + (A z ) 2 = 2 




Thus, in one experiment the solution will look like as Fraunhofer distribution. 
Figure 7 presents the surface (A) 2 Z=Q of the solutions (|3Tj) - (j3^j) . 

We can recognize the vortex structures of these solutions only if we see 
the field distribution and the intensity of one of the component, for example 
(A x ). On Figure 8 we present the surface of one of the components (A x \ z= q), 
rewritten in Cartesian coordinates: 



(40) 



^x sin ( Jx 2 + y 2 + z' 2 

W = ^ 2 + 2 " 2 

x z + y z + z z 

It is clearly seen, that the surface of our three dimensional solutions possess 
linear phase dislocations in space, which characterizes the vortex structures 



The intensity of this component (A 2 ), which is one I = l,m = 1 dis- 
tribution, is presented in Figure 9. This picture is possible to see in one 
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Figure 8: Surface of the one of the component (A Xt \ z=0 ) of the averaged 
with one time period vortex solutions (|40p. It is clearly seen linear phase 
dislocation for these solutions. 



experiment, when an optical vortex of kind (|36 p -([38| ) is passed trough the 
linear polarizer. Then the internal structure of these vortex solutions com- 
posed by spherical harmonics with / = l,m = ±1,0 distribution, will be 
seen. Our solutions (pjjl -lpIH j) admits total angular momentum I = 1 and 
line phase singularities for the different components. That is why, in spite of 
the fact that the intensity field is total symmetric, we recognize the internal 
vortex structures of these solutions. 

6 Conclusion 

In conclusion, we investigated the application of the nonlinear amplitude 
equation (0), in different spectral zones, starting from the optical region, 
and extending up to the R6 region of dielectrics. Different generalizations 
of the paraxial equation are obtained, depending from the spectral zones, 
the material constants, and the initial shape of the pulses (long pulses or 
optical bullets). We find large spectral region near the plasma frequency in 
cold plasma for fs and few picosecond optical pulses, where the temporal 
effects of second order in the amplitude equations (JTJ) can be neglected. The 
same result is obtained also for spectral regions on the 'wings' of the electron 
resonances in metal vapors. In the last case the conditions to obtain such 
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Figure 9: The intensity profile of one of the components (A Xj \ z=Q ) 2 , which 
corresponds to I = l,m — 1 distribution. 



effects depend strongly on the number density of the gases. In these spectral 
zones the propagation of optical pulses is governed by the amplitude vector 
nonlinear Schrodinger equation (J25J) and, as it was predicted in [TH], stable 3D 
optical vortices can be observed. They look like as Fraunhofer distribution 
in the space. After passing trough the linear polarizer the vortex structures 
will be seen. The vector theory presented in this paper can be extended also 
in respect to different kind of nonlinearities and polarizations of the optical 
field. 

7 Appendix 1. Nonlinear polarization of two 
and three components of the same frequency 

The electric field associated with linear or elliptically polarized optical wave 
can be written in the form 

E(x, y, z,t) = — {{xA x + yA y ) exp(iu t) - c.c.) , (41) 

where A x and A y are the complex amplitudes of the polarization components 
of a wave with the frequency u> . 
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7.1 Linearly polarized components 

Let investigate the polarization dynamics of two initially linearly polarized 
components of the electrical field. Then, there is no initial phase different 
between the components, and the complex amplitudes can be expressed as 
a product of two real amplitude functions with equal phases: $l(A x ) exp(i0), 
3? exp(i<ft). Calculating the relation between the linear polarized compo- 
nents in the same manner, as it was performed for plane wave j2Zj, we obtain 
the same kind relation for the components of the electrical field: 

Ey = $l{A y ) 

where $t(A y ) and dt(A x ) are real functions, not real constants as in the case 
of plane wave. We investigate the propagation dynamics of optical pulses 
with arbitrary localized amplitudes. The soliton case appear as partial case, 
when equal amount of nonlinear and linear phase shift are compensated, the 
initial phase different do not changed, and the linear polarization keeps on. 
For optical pulses with arbitrary amplitudes, as the media is isotropic, the 
diffraction and the dispersion add equal phases to the two components and 
they can be written generally in the form: 

A x = &(A X ) exp(i(K diff z + n disp t) + 0) = $l(A x ) exp(i$(z, t)), (43) 

A y = &(A y ) exp(t(K diff z + n disp t) + 0) = to(Ay) exp(t<S>(z, t)), (44) 

where K di ff is the phase shift due to diffraction, VL disp is the phase shift due 
to dispersion, $(z, t) is the equal phase components for the both amplitudes, 
and dt(Ay), $t(A y ) are the real part of solutions for the amplitudes on distance 
z in isotropic media. For arbitrary solutions of kind (J43|) . (|44|) the relation 
(|42|) keeps on for the real part of the amplitudes. It is following from the 
fact, that the diffraction and the dispersion effects do not change the initial 
linear polarization of the pulses (do not yield linear birefringence in isotropic 
materials) . The nonlinear polarization P n i in the case of isotropic medium 
is well known and can be written in the form |2*H1 l2*o] : 

P nl = A(u) (E -E*)E + (E ■ E) E*. (45) 
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Substituting the expression for the two linearly polarized component 
and (|44|) in (|45|) . right way we obtain for the nonlinear polarization of the 
different components [2*Ej : 



^=U^) + ^)(I^I 2 + I^I 2 )^, (JO) 



We point here, that the expression is valid not only for nonlinear nonresonant 
optical response, when A = B = const, but also for nonlinearity near to some 
of the resonances, when A(u) 7^ B(u). In the same way we can obtain, that 
the components of the nonlinear polarization for one linear polarized vector 
in 3D space take forms: 



(48) 
(49) 

Kx = \A{u>) + (|A,f + \A y \> + \A z f) A z . (50) 

7.2 Elliptically polarized components 

In the case of elliptically polarized light the relation ()42|) take form: 

E y _ 3t(Ay) 






2 + 




2 + \A z \ 2 ) 


) (1^1 


2 + 


\Ay\ 


2 + \A z \ 2 ) 


) (\A* 


r+ 


\Ay 


\ 2 + \A z \ 2 ) 



E x $1{A X ) 



exp(i(6 y -5 x )), (51) 



where 5 y — 5 X ^ 0, nn, tt/2, is the different of the phase functions for the 
components of the vector field. Since the diffraction and the dispersion add 
equal amount of phase shift to the different components in isotropic media, 
the phase different keeps on again. If we substitute E x and E y with such 
phase difference in the expression for nonlinear polarization (|45|h and take in 
mind that A = B = const for case of nonresonant electronic nonlinearities, 
the next expression for the components can be found 
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Pni = const 



A J 



2 -\A 
3 



, (l , A x + ~(A* x A y )A y exp(2i(5 x - 5 V )) 



(52) 



P%i = const 



A,, 



3 1-4x1 



A y + -(A* y A x )A x exp( 



-2i{8 x - 5 y )) 



■ (53) 



The last terms in equations (f5*2*J) and (f5*3*|) present degenerate four wave mix- 
ing process (u± = cl> 2 = oo ). Since 5 X — 5 y ^ 0, this leads to periodically 
exchange of energy between the elliptically polarized components. Thus, 
such one exchange of energy can be observed in isotropic material as gases, 
metal vapors and crystals. When the phase different vanish (5 X — 5 y = 0) 
in the case of linearly polarized components, the amplitude functions ad- 
mit again equal phases, and the nonlinear polarization ()52|) and ()53|) can be 
transformed to the nonlinear polarization of linearly polarized components 
f!46|) and ()47|). This is one generalization of the results obtained for para- 
metric solitons in media [SHj. When there is no phase different between 
the components, the parametric therm appear as usual cross-phase term in 
the dispersion relations, and can be take in mind to determinate the right 
amplitude constants. The polarization dynamics in materials which admit 
modal birefringece is more complicate. In single-mode fiber, where the lin- 
ear birefgingent effect appear |31| . additional polarization mode dispersion 
exists (221 • This lead to random birefringence in optical fibers and random 
phase different in the last parametric terms in (|52j) and (J53j) . Our numerical 
experiment show that the parametric effects with random phases are not ef- 
fective, and can be neglected. One transitions to the Manakov system with 
coefficients equal to one in front of the cross-terms j2Ij is possible again for 
birefringent materials at elliptically angle 35° between the main axes [3*3] . 
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